arXiv:1506.07325vl [math.PR] 24Jun2015 


Population processes sampled at random times 


Luisa Beghin* Enzo Orsingher' 


Abstract 

In this paper we study the iterated birth process of which we examine the first-passage time dis¬ 
tributions and the hitting probabilities. Furthermore, linear birth processes, linear and sublinear 
death processes at Poisson times are investigated. In particular, we study the hitting times in all 
cases and examine their long-range behavior. 

The time-changed population models considered here display upward (birth process) and down¬ 
ward jumps (death processes) of arbitrary size and, for this reason, can be adopted as adequate 
models in ecology, epidemics and finance situations, under stress conditions. 

Keywords and phrases: Yule-Furry process, linear and sublinear death processes, hitting 
times, extinction probabilities, first-passage times, Stirling numbers, Bell polynomials. 
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1 Introduction 

We study here different compositions of point processes, say K(t ) := L(H(t)), t > 0, where L and H 
are mutually independent and can possibly be birth, death or homogeneous Poisson processes. These 
compositions arise in the analysis of different population models, when, in some specific experimental 
circumstances, the time t must be replaced by a stochastic process H(t), t > 0. In our analysis we 
concentrate our attention on the case where the time increments are unit-valued and thus are properly 
represented by point processes, such as linear birth and Poisson processes. In these cases the composed 
process can be regarded as a randomly sampled population model. 

The simplest case of the composition of two independent homogeneous Poisson processes (with 
different rates) has been already studied in |1T]. More recently, the iterated Poisson process has been 
considered in [B], 

Our investigation concentrates here on the distribution of the first-passage times, i.e. 

Tfc := inf{t > 0 : K(t) = fc}, k > 0, 

which, in this context, substantially differ from the upcrossing and downcrossing times, i.e. 

Uk := inf{t > 0 : K(t) > k}, 14 := inf{t > 0 : K(t) < k}, k > 0. 

Indeed, the composed processes K jump over (resp. under) any level k with positive probability. 
The hitting probabilities Pr(T), < oo) display different behaviors. For example, for the iterated 
birth process K(t) = B\(B 2 {t)) (with B\,B 2 independent linear birth processes), the probabilities 
Pr (Tk < oo) attain their maximal value at k = 3. In the linear death process at Poisson distributed 
time the hitting probabilities display an oscillating behavior, which can be perceived when the initial 
size of the population is sufficiently large. 

Furthermore, we observe that in all the cases considered here the hitting probabilities Pr(Tfc < oo) 
depend only on the parameter of the outer process. 
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Moreover, the processes K display sample paths with upward or downward jumps of arbitrary size. 
Other models with the same feature, such as the space-fractional Poisson processes or, in general, 
time-changed Poisson processes with Bernstein subordinators, have been analyzed in |13| . Multiple 
jumps are also displayed by the generalized fractional birth processes studied in [T]. As in this 
last work, this property is reflected in the form of the equation governing the state probabilities 
Pk(t) := Pr{/\(f) = k}, k > 0, where the time derivative p' k [t ) is shown to depend, also here, on all 
Pk-j{t), for j = 0 


In particular, our analysis concerns the following specific cases: 

• Bi{B 2 {t )) 1 t > 0 

(section 3) 

• B(N(t )), t > 0 

(section 4) 

• D(N(t )), t > 0 

(section 5) 

• D(N(t )), t > 0 

(section 6), 


where B\ and B2 are independent linear birth (Yule-Furry) processes, N is the homogeneous Poisson 
process, D is the linear death process and D is the sublinear death process. The latter is characterized 
by the fact that the probability of a further death in [t, t + dt) depends on the number of deaths 
recorded up to time t (unlike the linear case where this probability depends on the number of surviving 
individuals). The sublinear birth process was introduced in [5], while the sublinear death process has 
been investigated in jl2j . in a fractional context. 

The processes considered here can be useful to model the population evolution sampled at Poisson 
times. The iterated birth process can be more appropriate for the cases where the time separating 
occurrences is rapidly decreasing. This structure can be applied, for example, in some experimental 
studies of diseases or epidemic diffusions. 

Time-changed birth processes of different forms have been studied in J7j and [E], with various 
applications to finance. Time-changed Poisson and birth processes arise also in the study of fractional 
point processes (see, for example, m , 0, 0 and 0 )- 


2 Notations 

List of the main symbols 


Linear birth 
Poisson 
Linear death 
Sublinear death 
Iterated birth 

Linear birth at Poisson times 
Linear death at Poisson times 
Sublinear death at Poisson times 


PROCESS 

B a (t) 

Nx(t) 

D^t) 

D„{t) 

Z[t) = B a (B x (t)) 
X{t) = B a (N\(t)) 
y(t) = D^N x (t)) 
y(t) = D^N x (t)) 


STATE PROBAB. 


pf (t), 0 < k < n 0 
pf (t), 0 < k < n 0 

Qk (*)> k > 1 
Qk(t), k > 1 
9^ (t), 0 < k < n 0 
q\ (t), 0 < k < n 0 


where no > 1 is the number of the initial components of the population and A, fi, a > 0. 


3 The iterated linear birth process 

Let us first consider the composition of a linear birth (Yule-Furry) process B a := B a (t),t > 0, with 
one initial progenitor and birth rate a > 0, with an independent process of the same kind, with 


2 


parameter A > 0, i.e. 


Z(t) := B a (B x (t)), t > 0. 


(3.1) 


As a consequence of the definition (3.1), the initial number of progenitors of the iterated linear birth 
process is random, since Z{ 0) = B a { 1) a.s. 

The probability mass function of this process can be written, for any fc > 1, as 

«k(t) := Pr{Z(t) = k}=EPv{B a (B x (t)) = k\B x (t)} 

OO 

= ^I p r{B a (j) = k}Pr{B x (t) = j} 

3=1 


We observe that 


= e 


-At 




— Xt — OL 


^V Qi (i- e -^) (i- e - At ) 

3 =1 
fc-1 

E 

1=0 


k-1 




l J 1 - e-“( i+1 )(l - e~ xt )' 


Q?(t) = 


P —At—a 


q 2 (t) = qf(t) 


1 - e-“(l - e~ xt ) 
1 — e~ a 


1 — e 2a (l — e xt ) 


= Ee~ aBx{t) 


We prove now that the probabilities k > 1, t > 0 decrease as k increases, since we have that 


Qki^-Qk-lit) = e 


— Xt — OL 


l 

k-1 


k — 2 


( -e~ a Y 


— a.\k — 1 


(-e- Q ) 


k—2 

^VZ-i;i-e-“(*+D(l- e - At ) 1 1 - e-“ fc (l - e~ At ) 


r^ — Xt—C 


i- e ~ a y 


/ /c — 2 

^ v - 1 y 1 - e-^+^l - e" A *) 


k—2 


= —e 


— At—2a 


^(l-e- At ) m e" 2am ^ 


ra—0 


Z=0 


fe - 2 


(-l)'e 


i„ — aZ(m+l) 


= - e - xt ~ 2a (1 - e _At ) m e _2am (l - e - a ( m +i)^k -2 < Qj 


m=0 


for any k > 2. 

The sample paths of the iterated birth process display upward jumps of size larger or equal to one. 
For this reason the analysis of the first-passage time through an arbitrary level k > 2, i.e. 

Tf = inf {t > 0 : Z(t) = k}, 

is of a certain importance. Moreover, we shall prove that Pr{Tj? < oo} is strictly less than one, 
as it happens for the iterated Poisson process (see |T2]) and thus any level k can be avoided with 
positive probability, because of ’’multiple jumps”. We present the explicit distribution of Tjf in the 
next theorem. 


Theorem 1 The distribution of Tj? reads 

OO 

Pr{Tjf G dt}/dt = Xe- M - ak ^2je~ aj { 1 - [{e a - e~ aj ) k ~ l - (1 - e~ a °) k - 1 ] , (3.2) 

j =i 

for any k > 2, 0 < t < oo. 
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Proof. By a conditioning argument we get 


Pr{T# G dt} 


EPr{rf edt\B a (i)} 

k-1 

^Pr {B a (B\(t)) = k - h,Ba(B\(t + dt)) = k} 

h=l 

oo k— 1 

A dt YjY^ Pr{ B a (B x (t)) = k - h,B a (B x (t) + 1) = k\ B x (t) = j}Pr{B x (t) = j} 

i=1 h= 1 

oo k— 1 

Adt^ j p r{S A (i) = j} Y Pr{ B a (j + 1) = k\ B a (j ) = k - h} Pr }B a (j) = k - h } 

i=i i 

oo k— 1 

Xe- xt ~ ak dtJ2je~ aj {l - e~ xt y~ 1 Y 
t=t h=i 


h ^ (e“ - l) fc (l - e-“ J ) fe - h - 1 , 


which coincides with (3.2). ■ 

We note that, in the special case k = 2, we get 


\ __ p-a\ 

Pr {T? G dt} = — -- —Ldt = Ae At (l - e~ a )[qf (t)} 2 dt 

[1 - e-“(l - e~ xt )} 

Remark 2 By integrating (3.3} ), we get that 

OO 

Pr{Tjf < oo} = e~ ak Y e ~ aj [( e “ - e^f- 1 - (1 - e^f- 1 ] , 

3 =1 


(3.3) 


which can be alternatively written as 

OO OO 

Pr{Tf < oo} = Y - e"“ fc ^ e"“ J (l - e"^)^ 1 


i=2 


o =i 


= (1 - e~ ak ) Y e~ aj (l - e-^y- 1 - e~ ak - a (l - e"“) 


—a\k — l 


3 = 2 
oo 


= (1 - e~ ak ) Y Pr{B a (j) = fc} - e-“ fc Pr{R a (l) = fc}. 

i =2 


(3.4) 


In order to analyze some special cases, we supply also the following finite sum form of the hitting 
probabilities given in (3.3): 


oo k— 1 

Pr{Tf<oo} = e-« k Ye~ aJ Y 

j—1 r —0 

k—2 


k-1 

r 


^_ e -ajy p a(k-r-l) _ (_ p -aj 


(—e~ aj ) 7 


= e 


2: Z /h — l\ p -a{r+ 1) , x 

“"‘£( r )(-ir -1) • 

r—0 ' 7 


(3.5) 


It is now easy to show that the following relationships hold, for k = 2,3,4: 

Pr {Tf < oo} < Pr{7lf < oo} < Pr{7)f < oo}, 
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since we have that 


Pr {Tf < 00 } 

Pi-{T? < 00 } 
Pr {Tf < 00 } 


e _2a 



e _2a 

l-e~ a ' 

1 + e _ “ 



1 + e~ a 


p -2a 

1, i 1 - e_ “ 

)(1 + e 3 “) 


1 + e~ a + e~ 2a 


The following figures (produced by the software R) describe the behavior of the probabilities \3.ffy , for 
different values of a. We remark that the scales are not the same in different figures. 


Figure 1: Hitting times Pr {Tf( < 00 } 
a=M4 a= 1 /2 
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a = 3 



2 10 20 * 30 40 50 


The figures above show that the probabilities Pr{Xjf < 00 }, for k > 3, are decreasing functions 
which, for 0 < a < 1, vanish more rapidly than for a > 1. 

The rate of decreasing is slow, for all a, since, from (|3.4|), we get that 


E Pr{Tf < oc} 
k=2 


since 


22e- 0 *'52(l-e- 0 *) k - 1 -'£e 

j =2 fc=2 j =1 


OO OO 

-a(j+l) 


^(e" Q - e -“0'+ 1 )) fe - 1 


J=2 3 =1 


—ot(j+l)—c 


k =2 

1 - e~ aj 


1 - e~ a + 


(1 - e"“ J ')[l - e~ a - e - a ^ +1 ) - e - a (^ +1 )- c 


E 

3 =2 


(1 - «"“) \ E 


. i=2 


1 - e _ “ + e- Q (J +1 ) 

{1 - e~ aj )[l + e-^+V] 
1 - e _a + e -“<J+i) 


e _3 “(l - e"“) 
1 - e _Q + e _2 “ 


D -3a 


= OO, 


v (l-e-^')[l + e- a ^' +1 )] 1 

inn -r—rrr— = -. 

j —>00 1 — e a + e ) 1 — e a 
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4 Linear birth process at Poisson times 

We now pass to the analysis of the process defined as 

X(t) := B a (N x (t)), t> 0, 

obtained as a composition of a linear birth process B a with an independent homogeneous Poisson 
process N x (t),t > 0, with rate A > 0. The process X can be regarded as a linear birth process 
randomly sampled at Poisson times. We first note that, for k > no, 


Qk,n 0 W : = Pr{X(t) = k\X(0) =n 0 } 


(4.1) 


3=0 

k—n o 


A: — 1 
k — no 


e ~ajn 0 _ e ~aj\ k n ° 

\ 1 j\ 


= E 


1=0 


k — no 

l 


(-l) / exp{At[l-e- a(no+I) ]}. 


In particular, for k = no, we get 


_ — \t(i- e-“"°) 


(4.2) 


(^) ^ 

and thus the waiting time for the appearance of the first offspring is exponentially distributed with 
parameter A(1 — e~ an °). It is easy to show that 


EX(t) =n 0 e A ‘ (e “" 1) . 

The variance of X can be obtained as follows 

VarX(t) = VarE[B a (N x {t))\N x (t)]+EVar[B a (N x (t))\N x (t)] 


(4.3) 


(4.4) 


= Var 


n 0 e 


*N x (t) 


■E 


n 0 e 2aNxW - n 0 e aNxW 


= no (no + l)e 


Xt(e 2a — 1) _ n 2 e \ i(e“-l) 


- n 0 e 


2At(e“-l) 


For the sake of simplicity, from now on, we assume no = 1 and denote for brevity q£i as q);. The 
factorial moments of X can be evaluated as follows: 


E {X(t) [X(t) - 1] ...\X{t) - r + 1]| X{t) = 1} 

oo 

= J2Hk-l)...(k-r+l)q£{t) 

k—r 

= [by@] 

1 (a ty 


(4.5) 


= r\e ~ xt ' 




7—0 1=0 ' / 


— rv'i\ r ~^ 1 


At g„ ai (A ty (l-e-^) 


r— 1 


= r\e ■- } y e 

3=0 


j\ [1 - (1 - e-<*i)] r+1 
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oo 




r— 1 


j—0 

oo 


r- 


^ (Xte ar y ^4 /r — 1 


i=o 


J! 


m—0 


m 


(-i) T 


r—1 .. ,\ 

= j(-l) m exp{-At(l-e a(r - m) )}. 

m=0 ' ' 


By using result (4.5) we can check formulas (4.3) and (4.41, by means of independent calculations. 

Our main interest here is to study the distribution of the first-passage time for an arbitrary level 
k , i.e. 

Tg = inf{f > 0 : X(t) = k}, 

for k > 2. In the next theorem we derive the following explicit distribution of T*. 

Theorem 3 For t > 0 and k > 2 we have that 

k ~ 2 fk — 1 \ 

Pr{T? e dt\X(0) = l}/dt = Ae" Qfc ^(-l)M \( e a(fe-i-b _ i)exp{-At[l-e-“ (i+1 )]} (4.6) 

1=0 ' ' 


Proof. We start by considering that 
Pr{T*edf|T(0) = l} 


fc-1 


= ^2Pr{B a (N x (t)) = k-h,B a (N\(t + dt)) = k} 

h=l 
k- 1 

= ^ Pr{5 Q (7V A (*)) = fc - fc, 5 a (7V A (t) + dN x (t)) = k} 

h= 1 

oo k— 1 

= MtJ2J2 Pr < B ^x(t)) = k — h, B a (N\(t) + 1) = k\ N\(t) = j} Pr {N x (t) = j} 

j=0 h=l 

°° (\+\j 

= A e~ xt dt a(j) = k- h,B a (j + 1) = k} 

j =o d ■ /j =1 


= ^ e" a(fe - ,l) (l - e"“) 

j =o fe=l ' ' 

k —1 /, k-h- 

= Ae" At -“ fc diW ■ 1 )(e“-lf £ 

fc=l ' ' z=o 

^ / j -i \ h 

= Ae- At -“ fc d<W '7 W-l)" £ 

h=l ' ' 1=0 

k —2 

= Ae- A *- afc dty](-l)'exp{Afe- a(Z+1) } ^ 

Z=0 h =1 

= Ae-“ fe dt g (* ; (-1)' exp{-At(l - c-^ 1 ))} “gf ( 

Z=0 ' ' h= 1 ' 


—a\h 


k h 1 (Xte-^+Vy 


l 


k — h — 1 


j=o 


r- 


1=0 

k-l-l 


l 


(-l) z exp{A te-“ (z+1) } 


A; — 1\ (k — h — 1 

l 


k-l-l 


'k-l-l 
h 


(e“ - l) h 
(e“-l)\ 


which coincides with (4.6). 
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Remark 4 In the particular case k = 2, formula (f-6) reduces to 


Pr {T? e dt | A’(O) = l}/rft = Ae"“(l - e"“) e - At(1 - e 

= e~ a Pr{ZG dt}, 

where Z is an exponential r.v. with parameter A(1 — e ~ a ) and this clearly shows that for the first- 
passage time through the level k = 2, the following result holds 


Pr {Tf < oo| X(0) = 1} = e~ a < 1. 
Remark 5 For any k > 2, by integrating formula we get 

Pr{T* <oo|*(0) = l} 


= e 


k — 2 

-“‘EH) 1 


/—0 
k—2 


k- 1\ e a ( fc -i-0 - 1 


l ) 1 — e _ “( z + 1 ) 


^—otk 




z=o 

oo 


fc - 1 
l 

k—2 


= e 


— a.k 


c E e_Qro E(-1)' 


a(fc-l-I) _ ^ ' g—a(/+l)m 

ra =0 

/fc - i 


m=0 


Z=0 


i 


^Oi{k-l)—Oil(rn-\-l) _ aZraj 


- E 


, —a(m+l) 


(1 


_ p -a(m+l)^fe-l _ ^lg—a(m+l)(fc—1) 


m—0 


— Oik 


E 


m—0 


oo 

ra =0 


a ^2 _ g-amjfc-1 _ ^_ijfc-lg-am(fc-l) 

9 —as(m+l)^ _ g—a(m+l)^/c—1 _ ^-alCg-am^ _ g-amjfc-l 

oo 

= (1 - e~ ak ) E e”“ m (l - 

m=l 

oo 

= (1 - e _Qfe ) E Pr { B a(m) = fc}. 


m=l 


(4.7) 


(4.8) 


^4s m the case of the iterated birth process, the hitting probabilities (4-8) are not affected by the 
parameter A of the inner process (N\, in this case). Moreover, we remark that the relationship given 
in the last line is equal to the analogous one presented in \3.4]) (except for the additional term in ( 3 ./)[ ), 
which is due to the different starting point of the birth process with respect to the Poisson one): 

Pr {T* < oo| *(0) = 1} - Pr{Tf < oo} = Pr{R a (l) = fc}, 

for any k > 2. Thus, for the iterated linear birth process, the probability of reaching any level k in a 
finite time is strictly smaller than the corresponding probability for B a (N\(t)),t > 0. This difference 
decreases monotonically with fc. 


Remark 6 For k = 2 the first line of formula (4-8) reduces to ( -}•?). In the special case fc = 3 we 












have that 


Pr {Tf < oo X{B) = 1} 


-3a 


e“ - 1 
1 - e~ a 

„- a l + e~ 2a 

1 + e~ a 

nX 


e a + 1 - 


1 + e~ 


1 I P, ^ 

= Pr{ T£ < oo | X(0) = 1} - - — < Pr{T* < oo | *(0) = 1} < 1. 

We show now that, for any k > 2 the distribution of the first-passage time for the level k 
ically decreasing: we present the following heuristic argument 

Pr {T? < oo| *(0) = 1} - Pr {T^ < oo| X(0) = 1} 

OO OO 

= (1 - e~ ak ) Y e _am (l - e - am ) fe - 1 - (1 - e-"^- 1 )) Y e -“ m (l - e" am ) 


k-: 


m =0 
/»+oo 


m—0 


~ {l~e~ ak ) g -ax(i _ e ~ ax y~ x dx — (1 — e 


0 —a(fc-l) 


/O 


/■+0O 

) / e" aa; (l-e- 

J o 


1 — e 


—otk 


1 _ g-a(fc-l) 


ak a(k — 1) 

l-e" a n + .-. + e - "^ 1 ) 1 + ...+ e"“( fe " 2 ) 

T 


< 


a 

(l~e~ a ) 

a 

a 


l + ... + e" a(fc - 2) 

e ~a(k— 2) 


Ac — 1 
1 1 




k k — 1 

g-a(fc-l) 


0 -a(fe-l) 


< 0 , 


for k > 2. 

The plots of Fig.2 confirm the decreasing structure of the probabilities Pr{T* < ool 


Remark 7 From formula (4-8) we derive the following equality 

->x 

V. hA 

1 — e -Q:fe 


at , = Pr {T* < oo | *(0) = 1} Prjr*, < oo} *(0) = 1} 
9k 9k — 1 ’ -] „ — rvk ^ g— a(k— 1) 


= - ^ e" 2am (l - e" am ) fe - 2 < 0. 

m= 1 

It is clear from 0 that, for large values of k, we get 

Pr {T? < oo| T-(O) = 1} « P^T*, < oo| *(0) = 1} 

and thus the hitting probabilities slowly change with k as Fig. 2 confirms. 

By applying formula ( 4-9) we can prove also the following relationship 

Pr {T* < oo| *(0) = 1} Pr{T* 2 < oo| *(0) = 1} 


1 — e 


— ak 


1 — e -a(k-2) 


= - Y e~ 2am (l - e - am ) k ~ 3 (2 - e~ am ) 


m—1 


= 2 


'Pi'{T^_ 1 < oo| x(o) = 1 } 

Pr{Tf_ 2 <oo|*(0) = l}' 

1 — g-a(fc-l) 

1 — e -at(k- 2) 


E 

m—l 


e -3am( 1 _ e - 


is monoton- 


) k ~ 2 dx 


X(0) = i}. 


(4.9) 


(4.10) 


amjfe-3 
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Figure 2: Hitting times Pr{T* < 00 } 
oc=1/4 a = 1/2 



a=2 a=3 



Formula (4-10) shows that, for large values of k, 


9k 9k -2 — 2(9k -1 — 9k- 2 )- 


AT 




5 Death processes at Poisson times 

We here consider linear and sublinear death processes at Poisson times. The linear death process, 
with initial population of n 0 individuals and death rate p > 0, has a binomial distribution, i.e. 


“ Pr{r> M (i) = fc| £> m ( 0) = n 0 } = (^e^l - 0 < fc < n 0 . 


(5.1) 


The probabilities (5.1l satisfy the difference-differential equations 


(5.2) 


ftPk (*) = (t) + p{k + l)Pk+i(t) for 0 < k < n 0 — 1 

f t Pn 0 it) = Pnoi*) for k = n o 

withpf (0) = l{ k =n 0 }- 

Moreover, we denote the sublinear death process, with death rate p > 0, as D(t), t > 0 and its 
distribution is given by 


Pk (*) : = Pr {D^(t) = k £> m ( 0) = n 0 } = 


—o — k 


(1 




1 < k < no 


k = 0. 


(5.3) 


The probabilities pj: satisfy the following difference-differential equations 

f TtP° 0) = ~K n o ~k + 1 )Pk(t) + p{n 0 - k)p% +1 (t) 1 < k < n 0 

\ sPo’W = 9- n oPi (*) ^ = 0 


with pj? (0) = lr fc=no j.. In order to catch the probabilistic mechanism underlying (5.3) we should 
consider that 

Pr{Dfj(t,t + dt) — Dfj,(t) = — 1 D^ft) = k} = p(no — k + l)dt + o(dt), 1 < k < no- 
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We note that in the linear death process the probability of cancellation of an individual in [t, t + dt) 
is proportional to the number of existing components at time t. On the other hand, in the sublinear 
case, this probability is proportional to the number of deaths recorded up to time t. This property 
makes the sublinear process more suitable for describing epidemics and the diffusion of rumors. 

The expected value of D(t), t > 0 can be evaluated as follows 


E D(t) = e-^^(n 0 -fc)( 1 - e^ 4 ) 

0 —^t\ n o+l 


(5.4) 


fc=o 


n 0 

n 0 


1- (l-e 


_ 1(1 — e -^)l 

/i ' dt 


n 0 




.k—0 


1 - (1 - e"^ 


„- ut \n 0 +1 _ e -»tyo+l _ ( nQ + j) ^ 

no + (1 - e -^) n0+1 + e>lt (! - e"^)” 0+2 + 1 - e"* 


= n 0 + l-e^[l-(l-e-^)" o+1 ], n 0 > 1. 

It is easy to check that lim t _ s . +oc , E D(t) = 0, as expected. For the linear death process, the hitting 
time of the fc-th event, i.e. 

T k = inf { s > 0 : D^(s) = fc} 

has the following distribution 

P t{t^ g dt} = “ e -^)"°~ (fc+1 V(& + 1 )dt, t > 0, 0 < k < 

Thus we get 

Pr {T? < oo} = 1, 

for any k. On the other hand, in the sublinear case, the distribution of 


n 0 - 


Tj: = inf{s > 0 : D^{s) = k} 


reads 


Pr {Tj? G dt) = e _#rt (l - e~ l * t ) ,l0 ^ k+1) ufa - k)dt, t > 0, 0 < k < n 0 . 
Moreover we get 


r+oo 


Pr{T^ < oo} = / Pr{T^ G dt}dt 


= (no - k) Y 


riQ—k— 1 

- E 

r—0 


r =0 

no — k 
r + 1 
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no — k — 1 
r 

(-l) r = 1, 




e -tU{r+l) dt 


(5.5) 


(5.6) 


for any k. 


5.1 


Linear subordinated 


death process 


Let 

y(t) := D^N x (t)), t> 0, 


11 



where D^ is a linear death process with parameter /.i and N x is an independent homogeneous Poisson 
process with parameter A. Then, for 0 < k < no, the probability mass function reads 


Qk(t) := Pr{ 7(7 = k\ 7(0) = no} 


(5.7) 


1=0 


Pr{ D^l) = k\ D„(0) = n 0 } Pr{JV A (7 = 1} 

riQ — k 

E 

3=0 


(T) "E ("7 *)(-!)'■ «e{-"(i - ^ tt+l) )}. 


Since 


<&(*) = e" At(1 " e 


3 ) 


the waiting time of the first death of 7(7 is exponentially distributed with parameter A(1 — e _Am °) 
and has the same form of the corresponding probability of the subordinated linear birth process, which 
was given in formula (4.2). The probability generating function of 7(7 can be written as follows: 


= E 


ra—0 


1=0 

n 0 

m 


G y (u,t ) = e -At£M [ i_ e -^(i_ u )] 


(5.8) 


(-l) m (l - u) m exp{—At(l - 


The mean value and variance are respectively given by 

E7(7 = noe^*-*- 1 ) 


(5.9) 


and 


Vary (t) 


Var E [D^Nxit))] N x (t)} + E Var [D„(N x (t))\ N x (t)} 


Var 


n 0 e 




-E 


n 0 e 




( 1 - 


0 -p,N x (t) 


no (no - l)e xt & - n 2 e 2Xt ^~^ 


n 0 e 




(5.10) 


We are interested now in the differential equation satisfied by the distribution of the process 
7(77 > 0. As a preliminary result we prove the following lemma. 


Lemma 8 The probability generating function given in 
lem: 


satisfies the following initial-value prob- 


gG y (u,t) = - A 
G y (u,0) = u n ° 


1 _ e (l-«)(l-e 


G y (u, t) 


t,U> 0, 


(5.11) 


where e adu := J2T=o “¥7 ^ ie shift operator. 


Proof. By considering the equation governing the state probabilities of the Poisson process, we can 
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write 


A OO A 

—G y (u,t) = MNx(t) = l}G D (u,l) 


dt v ^ dt 

1=0 


= -\J2MNx(t) = l}G D (u,l) + \J2Pr{ N x(t) = l ~ 1 }G D (u,l) 
1=0 1=1 
oo 

= -XG y (u,t) + \J2MNx(t) =r}G D {u,r+l) 

r =0 
oo 

= —XG y (u,t) + A^Pr {N\(t) = r}e dr G D (u,r), 


r =0 


where G D is the probability generating function of D. Since, for any t > 0, 
^ G D (u,t) = n(l - u)-^-G D (u,t) 






G^fu.t) 


= -^ 2 (1 - +/i 2 (l - u) 2 -^G D {u,t) 


^G D (u,t ) = At 3 (l- M )|-G D (n,t)-3^ 3 (l-n) 2 |^G D ( U ,t) + 


03 

+/r 3 (l — u) 3 -j^G D (u, t ) 


and, analogously 
<9 4 


r) r) 2 

^ G D (u,t ) = -/r 4 (l - u) — G D {u,t) + 7//(l - u) 2 — G D (u,t) + 

r) 3 f> 4 

-6/(1 - M ) 3 ^G D (n,t) +/(l - u) A -^G D {u,t), 


we can write 




1=0 ^ J ' 


where {^} are the Stirling numbers of the second kind (see [H]). Thus we get (by considering that 

{*} = 0) 


8[ 

dt 


OO OO j J / . N Ajf 1 

G y {u , t) = -XG y {u, t) + A P*{Nx{t) = r} ^ ^ ^ l J _ U-1) J (1 - u) 3 ~ l —— T G D (u, r) 

r=0 7=0 1=0 3 


= -A G y (u,t) + A^Pr {N x (t) = r}£ " uJ^G^u.r) 


r—0 


7! 

j '=0 ^ ;=o 




= -AG y (w,f) + A 

o=o J ' 


(U - 1} ^ 


G y {u,t), 
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where we have applied the Dobinski formula for the Bell polynomial B ni i.e. 


B n (x) = 


k =0 


Equation (5.11| follows by applying the well-known exponential generating function of the Bell’s 
polynomial (see e.g. [H]), i.e. 

E = «p{(e* - 1),). 


n —0 


n\ 


Theorem 9 The probability mass function ofy(t)),t > 0 satisfies the following equation 

no — k 


d 

dt 


with the initial condition 


l 7,\ 

(*) = Ae-" fc - e" M ) r - Ag^(t), 0 < k < n Q , (5.12) 

r=0 ' ' 

1 = <Z ■ 


Proof. In order to derive the differential equation satisfied by (5.7) we rewrite the first line in (5.11) 
as follows: 


§- t G y (u,t) = ^E 11 e —^ d’ u G y (n,t)-\ G y (u,t) 

= A E atM E (1 ~ “ >,( ' l ~ — i» k - AC y («,» 


k —0 j =0 

no oo 


j- 


= —— k(k- j + l)u k - XG y (u,t) 


k =0 j =0 

no k 


1=0 


= A E^wE(i- e_M ) j y ufc_J E (,)(-«)' - *G y M 


k =0 j—0 

= [l 1 = k — j + l\ 

no k /, \ k 

k 


1=0 


k =0 j—0 

n o k 


= E -\G y (u,t) 


V=k—j 


k =0 1=0 ' ' j=k—l ' J ' 


k=0 1=0 

n 0 n 0 


3=k- 

k 


= aE“'E(,K<*) E( 1 - e "") J ( i _ t+;j )(-i) , -‘ +, -Ao 5 ’(«,t) 


,i 

1=0 k=l N 

= [j = fc -1 + j’} 


j=k—l 
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/j \ l / -I \ 

= A E E (i) ^(t) £(!- e~ fi y' +k ~ l ( ; j (-1/ - AG*( U , t) 

Z=0 A;—/ ' ' j v =0 ^ ' 

n o / j \ 

= A E u 'E f , jflfc W(1 - e-n k ~ l e-^ 1 ~ XG y (u, t) 

1=0 k—l ' A 

= [fc = r + Z] 

^0 / j_ 7 \ 

= A E U 'E( T ) q? +l (t)(l-e-») r e-» l -\Gy(u,t). 

1=0 r =0 ' ' 


Remark 10 By comparing equation \5.12 ) with (5.2), we can see that, in the subordinated case the 
time-derivative of (t) depends on the probabilities q~ ( t), for any k < r < no, while, in the standard 
case, the derivative of p(( depends only on p^ +1 . This is due to the fact that the process y(t),t > 0 
performs downward jumps of arbitrary size, whose law can be derived from (5.1) and can be written 
as follows: 


Pr{;y(f + dt ) = k | y(t) = r + k} = 


1 — Xdt + Xdte ^ k , 
\dt( r + k )(l - e-v) r e~» k 


r = 0 
, r > 1. 


(5.13) 


The first line in \5.13 ) means that, during the infinitesimal interval [t,t + dt), either the Poisson 
process (the ’’time”) does not change or it changes and all the individuals survive during a time span 
of unit length. 


In our view the main result of this section is the density of the first passage-time through the level 
k , i.e. 

Tfy = inf{s > 0 : 3^(s) = k} k = 0,... n 0 , 
which is presented in the following theorem. 


Theorem 11 The probability density of reads, for any k = 0, ...,no — 1, 
Pr{T* G dt\y(0) = n 0 }/dt 


(5.14) 


= \e~ xt ~^ k 


/ \ o° 


(A t) j e~^ k 


j! 


1 — e 


—m(j+!) 


riQ—k 


-[i- 


— fij ] n ° ^ 


Proof. We start by considering that, for k = 0, ...,no — 1 

Pr{T£ edt\y(0)=n o } (5.15) 

n 0 

= ^ Pl '{ D n(N\(t)) = h,D^(N x (t + dt)) = k) 

h—k-\-l 

np 

= J2 P*{Dn(N X (t)) = h,D ll (N x (t)+dN x (t)) = k } 

h=k -\-1 

oo no 

= A dtY^ E Pt { D hU) = h , D nU + 1) = k}Pv{N x (t) = j} + 

j= 1 h=k +1 

+XdtPr{D ti (0) = n 0 ,D ll (l) = k}Pv{N x (t)=0} 


15 













= t 


3=1 


r- 


h—k-\-l 


,-lik 


(1 -e~^y 


-^h-k n o\ -Mh 


*(l-e“ w ) 


-W\no — h 


+\dte~ M ( U ° l 


e '"'"(I - e M ) 


— !i\n 0 — k 


= £ 

+Ad£e 


j=i 

„-Xt( n O\-^k l 


h=k-\-l 

-fi\n 0 -k 


no — k 
no — h 


(1-e-^) 


-!i\h-k -Mh 


(1 — e~ M ) 


-fj.j\n 0 -h 


e -(l -e“T 

A dte- xt -» k (^^j J2 {[e" w '(l - e"^) + 1 - e -w]"°“ fc - [l - e -w] n °- fe j + 


j=i 

+Adfe _A * f ) e-^ fc (l - e -^) n °- fc 


which coincides with (5.14). 


Remark 12 Uy integrating (5.14) we get that, for k = 0, ...,no — 1 


3=0 


Pr{T% <oc\y(0)=n 0 } = e-» k l n °)Y / e ~'* jk \\ l - e ~ fl{j+1) } n0 * - [l - e~^] n ° * L (5.16) 


We note that the probability \5.16 ) has the same structure of formula (36) in l Ilf , which is related to 
the iterated Poisson process, despite the fact that the outer process D has decreasing paths. 


Remark 13 An alternative form of (5.14), as a finite sum, can be obtained as follows 

( \ ^ f Jc\ 

U °J ^ [°~ J(— l) m_1 (l — e~^ m )exp[—Xt(l — e^ k+m ))\, 

' m—1 ' ' 


which, by integration, gives 


Pr{ 7? < oo| y(0) = »„} = e-<“ M £(-1)”-' (”° *) y 

A ' m =0 A / 

for any k > 0. We note that, for k = 0, the extinction probability is given by 

n 0 / \ no 


-i f no — k\ 1 — e 


_ g— fi(m+k) ’ 


rap / \ JIQ 

Pr{T 0 ^<oo|^(0) = n 0 } = -^(-ir( °) = -£(-l) r 

r—1 ^ ' r=0 


i = l. 


For k = no — 1, we have instead that 


Pr{T^ o _ 1 <oo|y(0)=n 0 } = e~^ n ° ^no- 


1 - e~» 

_ g— u n o 


— g-/i(no-l) 


n 0 


< e 


1 + e~v + e~ 2 n + ... + e -/d™o-i) 

-li(no-l) n 0 _ i 


n 0 e 


-^i(no-l) 


Unexpectedly enough Fig.3 (which is obtained here for p = 1/2) shows that the probabilities 
Pr{ < oo| 3^(0) = no} do not display a monotonic behavior for sufficiently large values of no- 
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Figure 3: Hitting times Pr{T^ < 00 } 
n 0 = 10 n 0 = 20 



5.2 Sublinear subordinated death process 


Let _ ~ 

y(t) :=D^N x (t)) 

where D „ is a sublinear death process with parameter /x and N\ is an independent Poisson process 
with parameter A. Then the probability mass function reads 


9fc.no W := = k 3^(0) = n 0 } 


= ^Pr{£> M (Z) = k D m ( 0) = n 0 }Pr{JV A (t) = 1} 

1=0 

riQ — k 

~~ xt E 

j=0 


= e 


n 0 — k 
3 


(—l)- 7 exp jAte | 


for 1 < k < no, while, for k = 0, is given by 


9o!n 0 W = e At E °)(-l) J exp{Ate w }. 
3 =0 V J 


(5.17) 


(5.18) 


Clearly q$ nQ is a decreasing function of the initial number of individuals, as can be directly checked, 

by considering (5.31. The expected value of y(t), t > 0 can be evaluated, by considering (5.4), as 
follows 


E y(t) = e- At ^{n 0 + l-e w '[l- (l - e~^) no+1 ]} 

3 =0 J ‘ 

= n 0 + 1 - e- At(1 - eA1) + f](l - e — Ki)no+i ( Xte *y . 

3 =0 J ' 

< no-( l-e-^ 1 - 6 ^). 
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In the case of the subordinated sublinear death process, we have no markovianity and thus we cannot 
evaluate explicitly the distribution of the hitting times Tjj? , 0 < k < no — 1. For this reason we consider 
here the instant of the first downcrossing of the level k, i.e. 


= inf{s > 0 : 3^(s) < k}, 0 < k < no — 1. 


For k > 1, we can write 

Pr {V?>t} = Pr {y(t)>k} 

= [by ([5l7t] 
n 0 n 0 -l 


r—0 

The inner sum can be treated as follows 

riQ—r 


rio (LQ — 1 / 1 \ 

EE tT 

l—k r—0 ' ' 

lo-k riQ-r , 

<r*— n 7 — h ' ' 


E 

l=k 


no - l 
r 


riQ—r—k 


= E 


1=0 


no — l — k 


= [by eq. (2.25) in [2]] 
f no — k + D 
r + 1 


so that we get 


riQ — fc+1 

Pr {V?>t}= ]T 

r—1 


_ 1 (no - k + 1 


exp {—At(l — e /ir )} 


which vanishes, for t —> oo, for any k > 1. For k = 0, by applying (5.18) we have, instead, the 
following extinction probability: 


~ n 0 

Pr{V(? <t}= Pr {y(t) = 0} = e- Af £ 


3=0 



(—l)- 7 exp {Ate w }, 


which, for t — > oo, is equal to 1. 
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